Pseudorapidity Distribution of Charged Particles in d + Au Collisions 

at = 200 GeV 
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The measured pseudorapidity distribution of primary charged particles in minimum-bias d + 
Au collisions at ^ys NN — 200 GeV is presented for the first time. This distribution falls off less 
rapidly in the gold direction as compared to the deuteron direction. The average value of the 
charged particle pseudorapidity density at midrapidity is (dN c h/d»7}| J) |<o.6 = 9.4 ± 0.7(syst) and the 
integrated primary charged particle multiplicity in the measured region is 82 ± 6(syst). Estimates of 
the total charged particle production, based on extrapolations outside the measured pseudorapidity 
region, are also presented. The pseudorapidity distribution, normalized to the number of participants 
in d + Au collisions, is compared to those of Au + Au and p + p systems at the same energy. 
The d + Au distribution is also compared to the predictions of the parton saturation model, as well 
as microscopic models. 

PACS numbers: 25.75.-q, 25.75.Dw 
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The pseudorapidity distribution of charged particles 
in d + Au collisions is important for understanding the 
evolution of the system created in more complicated Au+ 
Au collisions, and may provide a constraint on the initial 
state parton density. The energy and centrality depen- 
dence of pseudorapidity distributions in Au + Au colli- 
sions measured at RHIC [J, @, 13 are consistent with the 
approach based on the ideas of parton saturation 0, [j| 
and semi-classical QCD |(| . Measurements of d + Au col- 
lisions may be the ideal way to search for the onset 
of gluon saturation, since the system should be much 
simpler than the one studied in Au + Au collisions. The 
results of d+ Au collisions are therefore crucial for testing 
the saturation approach . 

In this Letter we present the first measurement of 
the minimum-bias pseudorapidity distribution of primary 
charged particles (dN c h/d7y) produced in collisions of 
deuterons with gold nuclei at a nucleon-nucleon center- 
of-mass energy, y/s NN , of 200 GeV. The pseudorapidity, 
rj, is defined as r\ = — ln[tan(#/2)], where 9 is the emission 
angle relative to the direction of the deuteron beam. 
The data were obtained with the PHOBOS detector 
at the Relativistic Heavy Ion Collider at Brookhaven 
National Laboratory. The data were collected using 



the multiplicity array covering | r\ |< 5.4. The array 
consisted of a barrel of silicon detectors surrounding the 
beam pipe in the central rapidity region ("Octagon"), 
and six forward silicon counters, three on each side of the 
interaction point ("Rings"). The multiplicity array used 
in d + Au collisions was the same as that for Au + Au col- 
lisions at y/s NN = 200 GeV |3|. The detector setup also 
included two sets of 16 scintillator counters, covering 
3 <| r\ |< 4.5. These counters were used in the primary 
event trigger and in the offline event selection. 

The results presented in this Letter were obtained by 
two analysis methods: a "hit-counting" method and an 
"analog" method. The details of the analysis proce- 
dure of the two methods leading to the measurement 
of dN c h/d?7 can be found in Ref. pj. The measured 
pseudorapidity distribution was corrected for particles 
which were absorbed or produced in the surrounding 
material and for feed-down products from weak decays 
of neutral strange particles. Uncertainties associated 
with these corrections, which we take as 20-50% of the 
corrections, dominate the systematic errors. 

Due to the low multiplicity in d + Au, a new algorithm 
for collision vertex reconstruction was developed, using 
the hit position and energy deposited in the Octagon 
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FIG. 1: Minimum-bias pseudorapidity distributions of pri- 
mary charged particles, dN c h/d?7, measured for d + Au 
a t V s nn — 200 GeV. a) dN c h/d7? distributions obtained by 
explicit integration over centrality bins for each centrality 
measure, b) Final dN c h/dr7 distribution obtained by averag- 
ing over the five minimum-bias distributions. The gray band 
corresponds to the systematic errors (90% C.L.). The curve 
corresponds to a triple Gaussian fit to d + Au. 



detector. Monte-Carlo studies using HIJING ^(j an d 
GEANT 3.21, as well as comparisons with the vertex 
information from reconstructed tracks, show a resolution 
of (Tvtx(z) = 1.4 cm for the most peripheral events and 
Cvtx(z) = 0.8 cm for the most central events, z being the 
distance along the beam axis. 

In order to determine the minimum-bias pseudorapid- 
ity distribution for this analysis, it is necessary to correct 
for the trigger and vertex finding efficiencies. Since these 
efficiencies are functions of multiplicity (and therefore 
centrality), the most model- independent approach is to 
measure the pseudorapidity distribution in narrow bins 
of centrality and then integrate over centrality to produce 



a minimum-bias result. The alternative approach of mul- 
tiplying by a global Monte-Carlo (MC) correction factor 
(reconstructed/true MC) would rely on our model to get 
the details of the shape and the centrality dependence 
right. To further reduce the model-dependence, we used 
five distinct measures for estimating the collision cen- 
trality based on the multiplicity (defined by the energy 
deposited in the silicon detectors) in different regions of 
pseudorapidity. The first of these centrality measures, 
Exot, used data from the full pseudorapidity coverage of 
PHOBOS, | r\ |< 5.4. The second measure, Eo c t, used 
data from the more central region | r\ |< 3. The next 
two measures, EAuDir and EdDir, used data from the 
gold direction (—5.4 < 77 < —0.5) and deuteron direction 
(0.5 < rj < 5.4), respectively. The final centrality mea- 
sure, Epjj ng , used data far from the mid-rapidity region 
from the silicon ring counters at 3.0 <| r\ |< 5.4. 

The centrality measures were calibrated using HI- 
JING 10] and GEANT 3.21 simulations. In the case 
of E-Tot, Eoct, EAuDir and EdDir, the distributions in 
data and MC were very similar once the triggering and 
event selection were applied to the MC and an additional 
scaling factor near unity was included. For instance, for 
ETot this scaling factor was 1.046. This allowed us to use 
the zero-bias "true" distribution in the MC to divide the 
data into equal 10% cross-section bins and to estimate 
the overall trigger and vertex-finding efficiency of ~ 83%. 
Based on studies of reconstructed and true pseudorapid- 
ity distributions from different Monte-Carlo simulations, 
we include a global 5% systematic uncertainty on the 
pseudorapidity distribution in the reported errors. In 
the case of ERi ng , the MC and data distributions were 
sufficiently different that the procedure had to be modi- 
fied. From the calculated efficiencies in each of the MC 
centrality bins, we could derive how many of the events 
in each "true" 10% bin would appear in the data. Cuts 
were then made on the Epu ng distribution that gave this 
efficiency-corrected number of events. 

Once the efficiency had been calculated and central- 
ity cuts made, we extracted the average number of 
participating nucleons (N part ) for each centrality bin 
using HIJING+GEANT. The measured minimum-bias 
pseudorapidity distribution and the estimated number 
of participants were obtained by averaging over all cen- 
tralizes for each centrality measure. This averaging was 
done as the minimum-bias distribution is less sensitive 
to centrality resolution and multiplicity bias effects than 
the individual centrality bins. The final distribution 
was obtained by averaging over the five minimum-bias 
distributions given by the different centrality measures. 
The total number of participant nucleons, {N part }, was 
estimated to be 8.1 ± 0.7(syst). To cross check the final 
distribution obtained by using the silicon centrality mea- 
sures, a further analysis was performed, only requiring at 
least one hit in one of the scintillator counter arrays. This 
analysis was consistent within the present systematic 
uncertainty. 

Figure shows the minimum-bias pseudorapidity 
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FIG. 2: dN c h/dri distributions per participant are shown for 
minimum-bias d + Au collisions at ^/s NJV = 200 GeV (solid 
points), for central Au + Au collisions (0-6%) jgj (triangles) 
and for p + p collisions from UA5 [l6j (open squares) at the 
same energy. The systematic errors are shown by gray bands. 
The uncertainties on {N pa rt) for d + Au and Au + Au have 
been added in quadrature to the gray bands. 



distributions of primary charged particles measured 
for d + Au collisions at y/s NN — 200 GeV obtained with 
the five centrality measures described above. The distri- 
butions agree to within 5% everywhere. Fig.^D shows the 
final result obtained by averaging with equal weights over 
the five distributions in Fig. The systematic errors 
(90% C.L.) are shown as a gray band. The statistical 
errors are negligible. Fig. ^> clearly shows that the 
cross section for particle production is largest in the gold 
direction (r\ < 0), and smallest in the deuteron direction 
(r\ > 0). By fitting the gold-side using a Gaussian in the 
region —3 < r\ < —0.5, the peak centroid was found to be 
at rj ~ — 1.9. This peak position is much closer to mid- 
rapidity than predicted by the parton saturation model 
Q, see Fig. 3a. ft may be interesting to note that a 
double Gaussian decomposition of the Au+Au distribu- 
tion at ^/s^ = 200 GeV Q for the — 6% most central 
collisions places the Gaussian centroids at rj ~ ±1.9, see 
Fig. 2. The significance of this observation is, however, 
not obvious, and may indeed be accidental. 

For d + Au, the measured average value of the 
charged particle pseudorapidity density at midrapidity 
is (dN c h/d^)| ?; |<o.6 = 9.4 ± 0.7(syst) and the integrated 
primary charged particle multiplicity in the measured 
region is Nj£, <s 4 = 82 ± 6(syst). Since dN ch /d?7 falls 
off to only about half of its maximum value at the 
most negative pseudorapidity of the detector acceptance, 
it is non-trivial to estimate the total charged particle 
multiplicity, N c h, for this system. However, using a 
triple Gaussian fit to the data we obtain the solid curve 
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FIG. 3: Comparison of the measured minimum-bias 
pseudorapidity distribution for d + Au collisions at 
^ys NN — 200 GeV to model predictions, a) Comparison to 
the parton saturation and RQMD models, b) Com- 
parison to the predictions of the HIJING [lOj model and the 
published AMPT model with and without final-state inter- 
actions (FSI) [l^|. c) Quantitative evaluation of the model 
predictions, expressed as the ratio of the model prediction to 
the data. The gray band corresponds to the systematic errors 
(90% C.L.) on the data. 



in Fig. ^d, which corresponds to a total number of 
produced charged particles of N{rjj t = 99 and the fit of 
the upper limit of the systematic errors (upper limit of 
the gray band) gives N c h = 110. We believe, however, 
that this represents an extreme upper limit for N c h since 
experimental data onp+A multiplicity distributions at 
lower energies 0, 0, 0], as well as various 
theoretical models discussed below, indicate that the fall- 
off at large negative pseudorapidities is more abrupt than 
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indicated by this fit to the data. A more reasonable 
extrapolation, guided by the expectations based on lower 
energy p+A data and predictions of the AMPT model, 
yields an estimate of N c h = 5 in the extrapolated region. 
Therefore, we estimate that the total charged particle 
multiplicity in this reaction is N c h = 87^Qy. It should be 
kept in mind, however, that even this estimate contains 
a significant contribution from the unmeasured region of 
pseudorapidity. 

Figure [21 shows the normalized minimum-bias 
dN c h/d?7 distribution for d + Au in comparison to 
the most central (0-6%) Au + Au collisions Q and 
non-single-diffractive p + p collisions from UA5 0] 
at the same energy, ^/s~^ = 200 GeV. To be 
consistent with our previous Au + Au pseudorapidity 
distribution publications, we have used a normalization 
factor of (Np ar t/2) which is unity for p + p. The 
distributions from different systems are normalized 
chiefly to enable a comparison of their shapes. The 
distribution seen in d+Au collisions at RHIC has 
features similar to those seen in p+A collisions at lower 
energies 0,0,0,0,0- Compared to p + p collisions, 
there is a significant increase in particle production in 
the gold fragmentation region and a reduction of particle 
production at positive pseudorapidity in the deuteron 
direction. The overall normalized production of particles 
is approximately the same as in proton-(anti-) proton 
collisions |lq |. 

Figure shows a comparison of the measured 
minimum-bias pseudorapidity distribution to the pub- 
lished predictions of the parton saturation model 
and the results from RQMD It is evident that 

these models are inconsistent with the data. The parton 
saturation model overestimates the height of the gold- 
side peak, underestimates its width, and predicts the 
peak at r\ ~ —3 rather than 77 = —1.9 as in the data. On 
the deuteron side, the model underestimates the charged 
particle production. RQMD drastically underpredicts 
the particle production and has surprisingly little pseu- 



dorapidity asymmetry. 

Figure presents the comparison of the mea- 
sured dNch/dry distribution to the predictions of HI- 
JING H3 and the published AMPT calculation [3. The 
HIJING calculation (dash-dotted curve) reproduces the 
deuteron side and the peak of the gold-side, but fails 
to reproduce the tail in the gold direction (77 < —2.5). 
The AMPT model uses HIJING for the initial stage of 
the collision, adding a " pop corn mechanism" for baryon- 
antibaryon production |l8l Il9| and final-state interac- 
tions (FSI). The solid curve is the default, while the 
dashed curve excludes FSI. We see that both final- 
state interactions and the popcorn mechanism appear to 
broaden the gold-side peak, leading to a moderate in- 
crease of the particle multiplicity in the region r\ < —3.5. 
The ratio of the model predictions to data is shown in 
Fig. 

In summary, the pseudorapidity distribution of 
charged particles produced in d + Au collisions 
at y/s NN = 200 GeV has been measured. The 
distribution is rather broad and peaked in the gold 
direction similar to observations in p+A collisions at 
lower collision energies. The measured pseudorapidity 
distribution is compared with the predictions of 
the parton saturation model as well as microscopic 
models. While AMPT predictions fall close to the data, 
the saturation model overestimates the asymmetry 
between particle production in the gold and deuteron 
hemispheres. 
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